In this notebook, we are using the maf_autopsy.tsv file
generated from the 01-preprocess-data.Rmd script. We will
also use the oncoprint goi list from OpenPedCan
oncoprint-goi-lists-OpenPedCan-gencode-v39.csv.
suppressPackageStartupMessages({
library(tidyverse)
library(ggrepel)
})
# Detect the ".git" folder. This will be in the project root directory.
# Use this as the root directory to ensure proper sourcing of functions
# no matter where this is called from.
root_dir <- rprojroot::find_root(rprojroot::has_dir(".git"))
scratch_dir <- file.path(root_dir, "scratch")
analysis_dir <- file.path(root_dir, "analyses", "tmb-vaf-longitudinal")
input_dir <- file.path(analysis_dir, "input")
# Input files
maf_autopsy_file <- file.path(scratch_dir, "maf_autopsy.tsv")
oncoprint_goi_file <- file.path(input_dir, "oncoprint-goi-lists-OpenPedCan-gencode-v39.csv")
# File path to plots directory
plots_dir <-
file.path(analysis_dir, "plots")
if (!dir.exists(plots_dir)) {
dir.create(plots_dir)
}
source(paste0(analysis_dir, "/util/function-create-corplot.R"))
source(paste0(root_dir, "/figures/theme.R"))
# Read maf_autopsy file generated from step `01-process-data.Rmd`
maf_autopsy <- readr::read_tsv(maf_autopsy_file, guess_max = 100000, show_col_types = FALSE)
# Read oncoprint_goi_file
# We will use the following list of genes and then we will create a label to select the ones for the corplots.
oncoprint_goi <- read.csv(oncoprint_goi_file, stringsAsFactor = FALSE) %>%
select(LGAT, Embryonal.tumor, HGAT, Other)
# let's create a list with the genes to be used for the corplots
oncoprint_goi <- data.frame(goi_list = unlist(oncoprint_goi))
# remove empty rows
oncoprint_goi <- oncoprint_goi[!apply(oncoprint_goi == "", 1, all), ]
oncoprint_goi <- data.frame(oncoprint_goi)
# We will define samples based on the "Kids_First_Participant_ID" column
samples <- unique(as.character(maf_autopsy$Kids_First_Participant_ID))
print(samples)
[1] "PT_23NZGSRJ" "PT_37B5JRP1" "PT_6N825561" "PT_CXT81GRM" "PT_DNAJYFZT" "PT_HFQNKP5X" "PT_KBFM551M" "PT_KTRJ8TFY" "PT_KZ56XHJT" "PT_MDWPRDBT"
[11] "PT_MNSEJCDM" "PT_PR4YBBH3"
for (i in seq_along(samples)) {
print(i)
maf_sub <- maf_autopsy %>%
filter(Kids_First_Participant_ID == samples[i])
timepoints_other_plots <- unique(maf_sub$timepoints_other)
timepoints_other_plots <- timepoints_other_plots[!timepoints_other_plots == "4"]
print(timepoints_other_plots)
timepoints_deceased_plots <- unique(maf_sub$timepoints_deceased)
timepoints_deceased_plots <- timepoints_deceased_plots[!timepoints_deceased_plots %in% c("1", "2", "3")]
print(timepoints_deceased_plots)
# Run corplot
for (t in seq_along(timepoints_deceased_plots)){
for (k in seq_along(timepoints_other_plots)){
fname <- paste0(plots_dir, "/", samples[i], "-", timepoints_other_plots[k], "-vs-", timepoints_deceased_plots[t], "-vaf-corplot.pdf")
print(fname)
p <- create_corplot(maf = maf_sub,
timepoints_other_plot = timepoints_other_plots[k],
timepoints_deceased_plot = timepoints_deceased_plots[t],
sid = samples[i])
pdf(file = fname, width = 10, height = 8)
print(p)
dev.off()
}
}
}
[1] 1
[1] "Diagnosis_BS_05S9WJW6"
character(0)
[1] 2
[1] "Diagnosis_BS_EJV0N3BX"
character(0)
[1] 3
[1] "Progressive_BS_E0S2Y0TS"
character(0)
[1] 4
[1] "Diagnosis_BS_HZNKSQ17"
character(0)
[1] 5
[1] "Recurrence_BS_EJ1H9PZY"
character(0)
[1] 6
[1] "Progressive_BS_FBJ516WW"
character(0)
[1] 7
[1] "Diagnosis_BS_9P4NDTKJ" "Diagnosis_BS_M0B42FPR" "Diagnosis_BS_M5FM63EB" "Diagnosis_BS_WYTDVC0Y"
[1] "Deceased_BS_J8EK6RNF"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KBFM551M-Diagnosis_BS_9P4NDTKJ-vs-Deceased_BS_J8EK6RNF-vaf-corplot.pdf"
Warning: `aes_string()` was deprecated in ggplot2 3.0.0.
Please use tidy evaluation idioms with `aes()`.
See also `vignette("ggplot2-in-packages")` for more information.
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KBFM551M-Diagnosis_BS_M0B42FPR-vs-Deceased_BS_J8EK6RNF-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KBFM551M-Diagnosis_BS_M5FM63EB-vs-Deceased_BS_J8EK6RNF-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KBFM551M-Diagnosis_BS_WYTDVC0Y-vs-Deceased_BS_J8EK6RNF-vaf-corplot.pdf"
[1] 8
[1] "Diagnosis_BS_3VKW5988" "Diagnosis_BS_402W79TS" "Diagnosis_BS_BQ81D2BP"
[1] "Deceased_BS_5968GBGT" "Deceased_BS_AF5D41PD" "Deceased_BS_EE73VE7V" "Deceased_BS_HYKV2TH9"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KTRJ8TFY-Diagnosis_BS_3VKW5988-vs-Deceased_BS_5968GBGT-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KTRJ8TFY-Diagnosis_BS_402W79TS-vs-Deceased_BS_5968GBGT-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KTRJ8TFY-Diagnosis_BS_BQ81D2BP-vs-Deceased_BS_5968GBGT-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KTRJ8TFY-Diagnosis_BS_3VKW5988-vs-Deceased_BS_AF5D41PD-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KTRJ8TFY-Diagnosis_BS_402W79TS-vs-Deceased_BS_AF5D41PD-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KTRJ8TFY-Diagnosis_BS_BQ81D2BP-vs-Deceased_BS_AF5D41PD-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KTRJ8TFY-Diagnosis_BS_3VKW5988-vs-Deceased_BS_EE73VE7V-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KTRJ8TFY-Diagnosis_BS_402W79TS-vs-Deceased_BS_EE73VE7V-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KTRJ8TFY-Diagnosis_BS_BQ81D2BP-vs-Deceased_BS_EE73VE7V-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KTRJ8TFY-Diagnosis_BS_3VKW5988-vs-Deceased_BS_HYKV2TH9-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KTRJ8TFY-Diagnosis_BS_402W79TS-vs-Deceased_BS_HYKV2TH9-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KTRJ8TFY-Diagnosis_BS_BQ81D2BP-vs-Deceased_BS_HYKV2TH9-vaf-corplot.pdf"
[1] 9
[1] "Progressive_BS_0ATJ22QA" "Progressive_BS_9DN4QR6E" "Diagnosis_BS_FWP8ZA4K" "Diagnosis_BS_H8NWA41N"
[1] "Deceased_BS_1Q524P3B" "Deceased_BS_22VCR7DF" "Deceased_BS_AK9BV52G" "Deceased_BS_D6STCMQS" "Deceased_BS_X5VN0FW0" "Deceased_BS_YHXMYDBN"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Progressive_BS_0ATJ22QA-vs-Deceased_BS_1Q524P3B-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Progressive_BS_9DN4QR6E-vs-Deceased_BS_1Q524P3B-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Diagnosis_BS_FWP8ZA4K-vs-Deceased_BS_1Q524P3B-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Diagnosis_BS_H8NWA41N-vs-Deceased_BS_1Q524P3B-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Progressive_BS_0ATJ22QA-vs-Deceased_BS_22VCR7DF-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Progressive_BS_9DN4QR6E-vs-Deceased_BS_22VCR7DF-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Diagnosis_BS_FWP8ZA4K-vs-Deceased_BS_22VCR7DF-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Diagnosis_BS_H8NWA41N-vs-Deceased_BS_22VCR7DF-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Progressive_BS_0ATJ22QA-vs-Deceased_BS_AK9BV52G-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Progressive_BS_9DN4QR6E-vs-Deceased_BS_AK9BV52G-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Diagnosis_BS_FWP8ZA4K-vs-Deceased_BS_AK9BV52G-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Diagnosis_BS_H8NWA41N-vs-Deceased_BS_AK9BV52G-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Progressive_BS_0ATJ22QA-vs-Deceased_BS_D6STCMQS-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Progressive_BS_9DN4QR6E-vs-Deceased_BS_D6STCMQS-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Diagnosis_BS_FWP8ZA4K-vs-Deceased_BS_D6STCMQS-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Diagnosis_BS_H8NWA41N-vs-Deceased_BS_D6STCMQS-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Progressive_BS_0ATJ22QA-vs-Deceased_BS_X5VN0FW0-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Progressive_BS_9DN4QR6E-vs-Deceased_BS_X5VN0FW0-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Diagnosis_BS_FWP8ZA4K-vs-Deceased_BS_X5VN0FW0-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Diagnosis_BS_H8NWA41N-vs-Deceased_BS_X5VN0FW0-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Progressive_BS_0ATJ22QA-vs-Deceased_BS_YHXMYDBN-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Progressive_BS_9DN4QR6E-vs-Deceased_BS_YHXMYDBN-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Diagnosis_BS_FWP8ZA4K-vs-Deceased_BS_YHXMYDBN-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_KZ56XHJT-Diagnosis_BS_H8NWA41N-vs-Deceased_BS_YHXMYDBN-vaf-corplot.pdf"
[1] 10
[1] "Recurrence_BS_00TRPEQX" "Recurrence_BS_BQWMQM2Y" "Progressive_BS_D368BNRD"
character(0)
[1] 11
[1] "Diagnosis_BS_CBMAWSAR" "Diagnosis_BS_ZSH09N84"
[1] "Deceased_BS_J8EH1N7V" "Deceased_BS_Y74XAFJX"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_MNSEJCDM-Diagnosis_BS_CBMAWSAR-vs-Deceased_BS_J8EH1N7V-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_MNSEJCDM-Diagnosis_BS_ZSH09N84-vs-Deceased_BS_J8EH1N7V-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_MNSEJCDM-Diagnosis_BS_CBMAWSAR-vs-Deceased_BS_Y74XAFJX-vaf-corplot.pdf"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/PT_MNSEJCDM-Diagnosis_BS_ZSH09N84-vs-Deceased_BS_Y74XAFJX-vaf-corplot.pdf"
[1] 12
character(0)
[1] "Deceased_BS_C2NH5FDT"